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In this paper, we show how a fixed point based boundary layer analysis can be used to under- 
stand phases and phase transitions in asymmetric simple exclusion processes (ASEPs) with open 
boundaries. In order to illustrate this method, we choose a two-species ASEP which has interesting 
phase transitions not seen in the one-species case. We also apply this method to the single species 
problem where the analysis is simple but nevertheless quite insightful. 



Asymmetric simple exclusion process (ASEP) [3 m 
its simplest form involves particles moving on a one- 
dimensional (ID) lattice with a bias in one specific direc- 
tion. Due to the mutual exclusion rule, the occupancy 
of a site by more than one particle is ruled out. The 
biased motion of particles is responsible for a particle 
current which is a typical indication of the system be- 
ing away from the thermodynamic equilibrium. With 
open boundaries, the system needs to be coupled to par- 
ticle reservoirs which maintain constant densities at the 
boundaries. All these non-equilibrium systems after suf- 
ficiently long time settle in a steady state where proper- 
ties like the average density of particles, current etc are 
independent of time. A striking feature of the open sys- 
tem is the boundary induced phase transition where the 
tuning parameters are the boundary densities 0. Var- 
ious bulk phases characterized through distinct current 
or particle-density profiles are usually represented in a 
phase diagram in the space of the boundary densities. 

Unlike equilibrium systems which carry no current, for 
ASEP, the information about the change in the bound- 
ary condition (BC) is mediated up to the bulk by the 
particle current. Therefore, the issue as how the BCs 
affect the boundary layers holds the key to understand 
the boundary induced phase transitions in ASEP. Recent 
studies 0, 01 have shown that the boundary layer analy- 
sis is a very general approach to probe different models 
of ASEP. In this approach, starting with the steady state 
hydrodynamic equation, a uniform approximation for the 
density valid both in the bulk and in the boundary re- 
gion is developed. For particle non-conserving ASEP, 
this method shows that the formation of a discontinuity 
in the bulk-density profile (or a localized shock) is due to 
a deconfinement of a boundary layer from the boundary. 
This deconfinement is shown to be associated with the 
divergence of the width of the boundary layer 0]. 

In this paper, we show how a fixed point description 
of the boundary layers and the phase-plane trajectories 
connecting the fixed points can be used to map the bulk 
phase diagram and identify the phase transitions. As 
an example, this method is applied to a particle number 
conserving ASEP that involves two species of particles 
In addition, we also consider the single species 
system, where the analysis is especially simple. 

In the two species model considered in this paper, par- 
ticles, denoted by A and B, move in opposite directions 



on a ID lattice of length L with N lattice points obey- 
ing the following hopping rules AO — > QA with rate 1, 
0-B BO with rate 1 and AB BA with rate 2. Here 
denotes an empty lattice site. While A type parti- 
cles are injected (withdrawn) at the left (right) boundary, 
the reverse happens for B type particles at these bound- 
aries. We assume that these injection and withdrawal 
rates effectively correspond to boundary reservoirs with 
fixed particle densities which are the BCs for our prob- 
lem. Two species models with specific boundary rates 
are known for spontaneous symmetry breaking among 
the two species even when their dynamics is completely 
symmetric. 0. While systems specified through bound- 
ary rates show spontaneous symmetry breaking, this is 
not seen in systems specified through boundary reservoir 
densities @. Our analysis is valid for arbitrary bound- 
ary densities. Domain-wall based studies with boundary 
densities [1| predict two first order transitions taking the 
system from a low to a high density jammed phase via 
an intermediate phase which is not seen in single species 
systems. In the steady state, the particle number conser- 
vation leads to a flat density profile in the bulk. Domain 
wall arguments show that unlike the high or low den- 
sity phase where, as in the single species case, the bulk 
densities are same as one of the boundary densities, the 
intermediate phase has bulk-density not directly related 
to the boundary densities. Thus, it is important under- 
stand what basic mechanism can be responsible for the 
existence of an intermediate phase and how the value of 
the bulk density in this phase is selected. The fixed point 
description presented here unveils these issues. 

The hydrodynamic formulation begins with the statis- 
tical averaging of the master equation that describes the 
time evolution of the probabilities of a given site being oc- 
cupied by either one of the two types of particles. A con- 
tinuum limit, as the lattice spacing, a — > and N oo 
with L ~ Na remaining finite, is, then, taken for the 
resulting equation. The differential equations, thus ob- 
tained, involve the average particle-densities p{x,t) and 
ri{x,t) for A and B type particles with x and t as con- 
tinuous variables representing position and time and the 
corresponding current-densities jp and as functions of 
p and 77. In the first order in e = a/ 2, the continuum 
equations are 



— -t- ^ = e— ((1 + 7])— ~ p—) 
dt dx dx dx dx 



(1) 
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The stationary state of the model on a ring is given by a 
product measure and the exact stationary fluxes, jp and 
jr,, are jp = p(l - p + 77) and jr, = -t]{1 -r] + p). 

The diffusion hke terms in equations (H]) and ^ are 
smaU in the e — > Hmit. However, these terms are crucial 
for generating uniform approximations of the solutions. 
The stationary bulk-densities {pi, and ?/;,) of constant val- 
ues, satisfy the differential equations trivially. These con- 
stant solutions cannot satisfy the BCs at both the ends. 
In order to satisfy the BCs, boundary layers are formed 
over narrow regions near one or both the boundaries. To 
focus on this narrow region, we need to scale the position 
variable {x — X{,)/e, where xq specifies the lo- 

cation of the boundary layer, and express the stationary 
equations in terms of this new variable. Here, xq ~ 0(1) 
for left(right) boundary layer which merges to the bulk 
particle-density profile asymptotically as a: — > -|-oo(— 00). 
A little algebraic manipulation and integrations of the 
stationary versions of equations ^ and ^ lead to 

= c + g(l — p), and 
[l+p) q{x)^d + p + q{c-pq), (3) 

where p = p+rj and q = p — rj with values lying within the 
range [0, 1] and [—1, 1], respectively. Here, prime denotes 
a derivative with respect to x and c and d are the inte- 
gration constants. The saturation of the boundary layers 
located, for example, near x ~ L, to bulk as i — > —00, 
requires c = — 9b(l — Pb) and d ~ q^ — pb, where pb and 
Qb are the bulk values of p and q, respectively. Since 
the bulk profile is flat, same c and d are applicable in 
case of boundary layers appearing near both the bound- 
aries. For particle non-conserving systems, where the 
bulk profile is x dependent, pb and qb are the values of 
p and q at the point where the boundary layer merges 
the bulk. Although this approach is applicable to parti- 
cle non-conserving systems, we restrict to the conserved 
case to avoid algebraic complications. Because of same 
c and d and the translational invariance, the same phase 
portrait is applicable to both the boundaries. 

To understand the origin of the three phases, it is use- 
ful to study the stability properties of the fixed points 
of equations in ([3]). A combined single first order differ- 
ential equation involving ^ leads to the solution, p{q), 
as a function of q, for a given BC at x = 1 or 0. These 
solutions are the phase-plane trajectories flowing to the 
appropriate fixed points. The boundary layers are these 
trajectories seen as functions of x. Since these fiows ap- 
proach the fixed points in i ±00, the bulk-densities 
correspond to one fixed point. Also, by the choice of the 
constants c and d, the equations are guaranteed to have a 
fixed point that corresponds to the bulk-densities. There 
exists three fixed points, (p* , q*), which are the solutions 
of the equations 

p* = {c + q*)/q* and {q*f-{d+l)q*-c^O. (4) 



Although the physically acceptable range of both c and d 
is [—1,1], there is only a narrow window in this c—d space 
(see the striped region in figure[T]), where the fixed points 
can represent a meaningful density profile (within the 
physical domain i.e. the values of p and q lie in the range 
[0, 1] and [—1,1] respectively and < p, 77 < 1 ). Outside 
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FIG. 1: Nature of the fixed points in the physically acceptable 
range of the c-d plane. On the bold solid line and the dashed 
line the discriminant of equation Q is zero. Fixed points are 
imaginary in the region outside these two lines. In the striped 
region enclosed by these lines and the lines d = c and d = — c, 
there are physically acceptable fixed points. In the rest of 
the space, not a single fixed point in the set of three has any 
physically acceptable value. 

this window, the fixed points are either imaginary, or do 
not have a single member that can represent the bulk. A 
stability analysis shows that out of the three fixed points, 
one is unstable, one is saddle and one is stable. Based on 
their locations on the p — q plane, we call the stable and 
unstable fixed points as right (R) or left (L). The saddle 
fixed point is always denoted as S. 
There exist the following possibilities: 

1. L is unstable and lies within the physical domain 
and R is stable with complex conjugate eigenvalues hav- 
ing negative real parts (spirally inward flow). R remains 
outside the physical domain and the bulk is, therefore, 
given by the L flxed point values. The bulk in this case 
is also the same as the left BC. The right BC is satisfied 
by a boundary layer which is a part of the trajectory on 
the p — q plane starting from L and reaching the stable 
fixed point R. R is outside the physical range of p so that 
any physical BC is satisfied by the right boundary layer. 
Trajectories corresponding to different BCs at a; = 1 are 
shown in figure (0\). 

2. R is stable and lies within the physical domain and L 
is unstable with complex conjugate eigenvalues (spirally 
outward flow) having positive real parts. By the above 
argument, here the bulk satisfies the right BC with a 
boundary layer at a; = satisfying the left BC. The L 
fixed point must lie outside the physical domain. The 
trajectories in this case are shown in figure (O?)- 

3. There exists another situation where the saddle 
fixed point and either one of R or L are in the physi- 
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FIG. 2: Flow trajectories to the stable fixed point (R) with 
the BCs as specified in figure. Arrows indicate the direction 
of increasing x. (A) c — .05 and d — —.65 with BCs spec- 
ified at a; = 1. (B) c = -.07 and d = -.41 with BCs at 
a; = 0. The rectangular boxes enclose portions of the trajec- 
tories with the BCs p(.21) = .95 for (A) and p(-.34) = .86 
for (B). The length of the box is fixed by the fixed point in 
the physical domain and the BC of the specific trajectory it 
encloses. The vertical width of the box is chosen in such a 
way that it encloses the relevant part of the trajectory of in- 
terest. These parts of the trajectories in A(B) are seen as the 
right (left) boundary layers in p-x and rj-x plane. 



cal domain but the latter is different from the BC. In 
this case, the boundary layers should be described by the 
separatrices such that the bulk densities can acquire the 
saddle fixed point values. BCs at both the boundaries 
are satisfied by boundary layers which arc parts of the 
separatrices joining L with S, and R with S. This phase 
where the bulk values of p and 77 arc given by the saddle 
fixed point values is the intermediate phase. 

A general inspection regarding the stability of the fixed 
points in the allowed range of the c — d plane confirms 
these rules as the only possibilities. It also indicates that, 
for c > (c < 0), the system remains either in L(R) or 
in the intermediate phase. In general, the intermediate 
phase appears when S is the only fixed point that lies in 
the domain specified by the BC. 

The three rules mentioned above completely decide the 
phases and the phase boundaries. Suppose, we are in the 
left phase where the bulk densities are same as the left 
BCs. This allows us to calculate c and d, and thereby. 



the explicit values of all the three fixed points. The fixed 
point corresponding to the bulk-density values, is an un- 
stable one and the boundary layer solutions approach 
these values as i — > —00 while approaching the stable 
fixed point, which lies beyond the physical region, as 
i — > 00. The boundary layers which arc part of the 
trajectories in figure appear as in figure ^ in the 
p-x or r/-x plane. The trajectories deviate because of the 




FIG. 3: Numerical solutions for the boundary layers in the L 
phase near x = 1 for c = .05 and d = —.65. The boundary 
layers of p and rj profile merge to the bulk values .2 and .7 
respectively. The solutions tends to saturate to the saddle 
fixed point given by {q,p) = ( — .153, .673) 

proximity of the saddle fixed point and this deviation ap- 
pears as the tendency of the boundary layer to saturate 
to an intermediate value. 

As the right BC is changed further keeping the left one 
unchanged, the solutions still appear like those of figure 
([2|) or ([3]) with same bulk-dcnsitics until a cusp appears 
in the trajectory. At this right BC, the trajectory is just 
the separatrix joining the unstable and the saddle fixed 
point on one side and the saddle and the stable fixed 
point on the other side. Since the fiow now meets the 
saddle fixed point, the bulk densities are given by the 
saddle fixed point values. This requires formation of a 
new boundary layer near a; = with the height p* — p* 
where p* and p* are the values of p at the unstable and 
the saddle fixed points, respectively. At this point, the 
system enters the intermediate phase. Since the values of 
the fixed points depend only on c and d or on the bulk- 
density in the L phase, the height of the newly formed 
boundary layer is only dependent on the bulk-density (or 
on c and d) of the L phase. The variation of the height of 
these new boundary layers with d for a given c is shown 
in figure 

The formation of a new boundary layer with a finite 
height and a discontinuous change in the value of the 
bulk density are indications of a first order phase transi- 
tion. This basic principle of the bulk densities acquiring 
the saddle fixed point values continues to hold good in- 
side the intermediate phase. In this phase, as the right 
BC is changed, the bulk values, ph and qh, change with 



4 



CD 




0.7-0.6-0.5-0.4-0.3 -0.2-0.1 



FIG. 4: Heights of the newly-formed boundary layers at a; = 
in the p and rj profiles are plotted with d for c = .05. These 
boundary layers are formed at the phase boundary between 
the L and the intermediate phase. 



boundary layers appearing at both the boundaries. The 
two boundary layers are the two separatrices joining the 
saddle fixed point to the stable and unstable fixed points. 
The bulk-densities or the values of c and d are such that 
the right and left BCs lie exactly on these two separatri- 
ces. 

Next, we apply this method to ASEP with only one 
kind of particles moving rightward on a ID lattice obey- 
ing mutual exclusion. In the steady state, there are three 
phases depending on the left and right boundary densi- 
ties a and 1 — /3, [l^l with the bulk particle-density pb 
being ph = 1 — /3, for a > /3, and /3 < 1/2, pb = a, for 
(3 > a and a < 1/2, and pb = 1/2, for a > 1/2 and 
/? > 1/2. The approach to these phases can be under- 
stood well through the motion of the domain walls [llj 
which separate possible steady states of the system. 

The equation that describes the boundary layer in this 
case is 



(l/2)p'(£) p(l - p) + c, 



(5) 



where the constant c is chosen as c = — Pb(l — Pb) to 
assure saturation of the boundary layer to the bulk as 
i — > oo or — oo. This equation has two fixed points p* = 
Pb, 1 — Pb- A stability analysis around these fixed points 
shows that (i) pb is stable for pb > 1/2, and (ii)l — pb is 
stable if pb < 1/2. 

In case of pb > 1/2, if a boundary layer has to saturate 
to Pb (the stable fixed point), its location should be at 
X = 0. Asi^— oo, this boundary layer saturates to the 
unstable fixed point 1 — pb- The right BC is satisfied if 
Pb = I ~ f3. Two conditions naturally follow from this, 
(i) 1-/3 > 1/2 or /3 < 1/2 and (ii) a > 1 - pt or a > /?, so 



that any boundary layer at x = satisfies the BC before 
saturating to the fixed point 1 — pb- Similarly, in case of 
Pb < 1/2, the boundary layer saturates to 1— pf, as 5: oo 
or ph as a; — > — cxD. A boundary layer of this nature should 
be located at x = 1 with the conditions {i)a < 1/2 and 
(ii) 1 — Pb > I — f3 or a < (3. As we see here, the location 
of the boundary layer is completely determined by the 
limit in which the bulk fixed point is approached by the 
solution. In both pf, > 1/2 and pb < 1/2 cases, while 
the boundary layers continue to be in the same position, 
their slopes change across the line a = 1 — /3. These 
boundary layers correspond to a different set of solutions 
of equation ([5]) and they approach the bulk fixed points 
at the same limits of x as before with opposite slopes. In 
general, the boundary layers saturating to two different 
densities are like localized domain walls of reference [ll| 
separating two possible steady states. For pb = 1/2, the 
two fixed points merge, and it is easy to see that the 
boundary layers from both sides of the system are able 
to approach the bulk density in a power law fashion as 
^ li-r either as i ^ cxd or as a; — oo. Thus, in this 

\x\ ' 

case, there are boundary layers at both the boundaries. 
From a stability-analysis like approach, it is easy to see 
that both the boundary layers have negative slopes. This 
implies that such a phase is possible when a > 1/2 and 
(3 > 1/2. These three possibilities cover the entire a-f3 
parameter space. 

In this paper, we show how a phase-plane analysis of 
the differential equations for the boundary layers can be 
used to predict boundary induced phase transitions in 
ASEP. This method is applied to single species as well as 
two species processes. We arrive at a general prediction 
that for flat bulk profiles, the bulk densities in different 
phases are given by various fixed point values of the dif- 
ferential equations describing the boundary layers. For 
two species case, we find out the bulk density of the inter- 
mediate phase which has no one species analogue. The 
fixed point analysis determines the nature of the transi- 
tions, the location of the boundary layers and how the 
boundary layers merge to the bulk densities. For non- 
constant bulk density, the fixed point of the boundary 
layer should match with the bulk-density value near their 
merging. In addition, we believe that this method can be 
applied to more complicated systems, such as those with 
more than two species of particles and the number of 
phase transitions can be predicted from the locations of 
the fixed points on the phase-plane. The generality of the 
method suggests its wider applicability to various other 
systems where boundary induced transitions can be seen. 
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